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The dramatic discovery with the Rossi X-Ray Timing Explorer satellite of remarkably coherent ~300-1200 Hz 
oscillations in the X-ray brightness of some sixteen neutron stars in low-mass binary systems has spurred theoret- 
ical modeling of these oscillations and investigation of their implications for the neutron stars and accretion flows 
in these systems. High-frequency oscillations are observed both during thermonuclear X-ray bursts and during 
intervals of accretion-powered emission and appear to be a characteristic feature of disk-accreting neutron stars 
with weak magnetic fields. In this review we focus on the high-frequency quasi-periodic oscillations (QPOs) seen in 
the accretion-powered emission. We first summarize the key properties of these kilohertz QPOs and then describe 
briefly the models that have been proposed to explain them. The existing evidence strongly favors beat-frequency 
models. We mention several of the difficulties encountered in applying the magnetospheric beat-frequency model 
to the kilohertz QPOs. The most fully developed and successful model is the sonic-point beat-frequency model. 
We describe the work on this model in some detail. We then discuss observations that could help to distinguish 
between models. We conclude by noting some of the ways in which study of the kilohertz QPOs may advance 
our understanding of dense matter and strong gravitational fields. 



1. INTRODUCTION 

It has long been expected (see, e.g., J7|,pl|) 
that important information about the intrinsic 
properties of neutron stars and stellar mass black 
holes, as well as about the physics of accre- 
tion onto them, could be extracted from their 
X-ray variability at frequencies comparable to 
the ~1-10 kHz dynamical frequencies near them. 
The Rossi X-ray Timing Explorer (RXTE) was 
specifically designed [|]]3£| to have the large area, 
microsecond time resolution, and high telemetry 
rate needed to probe this high-frequency regime. 

The value of access to these high frequencies 
was dramatically confirmed when, less than two 
months after the launch of RXTE, remarkably co- 
herent ~ 1000 Hz oscillations were discovered (see 
41 1) in the X-ray brightness of two neutron stars 
in low-mass X-ray binaries (LMXBs). These are 
the fastest astrophysical oscillations ever discov- 
ered. It is thought that they are being generated 
at or near the surfaces of these neutron stars. 



To date, brightness oscillations with frequen- 
cies ranging from ~300 Hz to more than 1200 Hz 
have been discovered in some sixteen neutron 
stars in LMXBs. High-frequency oscillations are 
observed both during thermonuclear X-ray bursts 
|pE-37 and during intervals of accretion-powered 
emission j4l],^]. They appear to be a character- 
istic feature of disk-accreting neutron stars with 
weak magnetic fields. 

In this review we focus on the high-frequency 
quasi-periodic oscillations (QPOs) seen in the 
accretion-powered emission. The weak-field neu- 
tron stars in which these kilohertz QPOs have 
been discovered were studied extensively with 
EX OS AT and Ginga, and a detailed physical pic- 
ture was developed based on their 2-20 keV X- 
ray spectra and 1-100 Hz X-ray variability (see 
Q|lf,H,[ni). In tllis Picture, LMXBs with weak- 
field neutron stars fall into two classes, called the 
"Z" and "atoll" sources after the shapes of the 
paths they trace, over time, in X-ray color-color 
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Atoll Source 




Figure 1. Side view of an atoll source (top) 
and a Z source (bottom) . The dark shading indi- 
cates the accretion disk. The light shading indi- 
cates the hot gas in the magnetosphere and the 
central corona that surrounds the neutron stars. 
The arrows in the lower panel indicate the cooler, 
approximately radial inflow that is thought to 
be present outside the central coronae of the Z 
sources. From [28]. 



diagrams g§. Modeling |o|-§2| of the X-ray 
spectra of these sources indicates that both types 
are surrounded by a central corona with a scatter- 
ing optical depth ~3-10, and that the Z sources 
also have a cooler radial inflow (see Fig. [j]) . 

The six known "Z" sources have mass accretion 
rates comparable to the Eddington critical rate 
Me and inferred surface magnetic fields ~ 10 9 - 
10 10 G. They display two different types of lower- 
frequency QPOs. The ~15-50 Hz "horizontal- 
branch oscillations" (HBOs) have frequen- 
cies that increase steeply with accretion rate 



and are thought to be produced by the mag- 
netospheric beat-frequency mechanism |Q,|lJ,|2|. 
The ~4-8 Hz "normal/flaring branch oscilla- 
tions" (N/FBOs) fH are thought to be caused 
by radiation-hydrodynamic oscillations in the ra- 
dial inflow §0Jl|. 

The — 15 known "atoll" sources are both less 
luminous and more weakly magnetic than the Z 
sources, with accretion rates —1-10% of Me and 
magnetic fields - 10 7 -5 x 10 9 G. No QPOs with 
frequencies < 100 Hz have so far been detected 
in any of the atoll sources, with the exception 
of Cir X-l, which has a -1-30 Hz QPO with a 
frequency that varies with its brightness (see Q). 

Considerable theoretical effort has been de- 
voted to understanding the mechanisms that pro- 
duce the kilohertz QPOs. In this review we first 
summarize the key properties of these QPOs and 
then describe briefly the models that have been 
proposed to explain them. The kilohertz QPOs 
commonly occur in pairs, and the existing evi- 
dence strongly favors beat-frequency models. We 
discuss application of the magnetospheric beat- 
frequency model to the kilohertz QPOs. The 
most fully developed and successful model is 
the so-called sonic-point beat-frequency model, in 
which the higher frequency in a QPO pair is the 
orbital frequency of gas at the inner edge of the 
Keplerian disk flow and the lower frequency is the 
difference between this frequency and the spin fre- 
quency of the neutron star. We outline the sonic- 
point model and describe some of the calculations 
that have been carried out to explore it. We con- 
clude by discussing observations that could help 
to distinguish between models and noting some 
of the ways in which the study of kilohertz QPOs 
may advance our understanding of neutron stars, 
dense matter, and strong gravitational fields. 

2. CHARACTERISTICS OF KHZ QPOs 

The kilohertz QPOs observed in the neutron- 
star LMXBs have strikingly similar properties 
fr4l|Ji2|| . These properties include: 

High frequencies. — The frequencies of the kilo- 
hertz QPOs range from -300 Hz to -1200 Hz. 

High maximum amplitudes. — The strongest 
kilohertz QPOs have fractional rms amplitudes 
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as high as ~15% of the total 2-60 keV countrate. 
The amplitudes of the kilohertz QPOs are much 
larger in the atoll sources than in the Z sources, 
where they are typically < 1%. 

High coherence. — QPOs with quality factors 
Q = i/Qpo/A^Qpo ~ 50-100 are common; in 
some sources, the kilohertz QPOs have Q values 
as high as 200. 

Frequency-M correlation. — In the Z sources, the 
frequencies of the kilohertz QPOs increase steeply 
with M. In many atoll sources, the frequencies of 
the kilohertz QPOs increase steeply with coun- 
trate (in at least some, the frequency-countrate 
relation shifted significantly between observations 
made just a few days apart |p|,p3||). Although the 
magnetic field strengths and accretion rates of the 
atoll and Z sources are thought to be quite dif- 
ferent, their frequency-countrate and frequency- 
accretion rate correlations are very similar. 

Similar frequency ranges. — The frequency 
ranges of the kilohertz QPOs are similar in the 
Z and atoll sources, even though these two types 
of sources are thought to have very different mag- 
netic field strengths and accretion rates. 

QPO pairs. — Kilohertz QPOs commonly occur 
in pairs and, except in Sco X-1 [[ll| (and possi- 
bly in 4U 1608-52 (§3[), maintain a constant fre- 
quency separation Av as they shift up and down 
in frequency by hundreds of Hertz. 

Relation to spin frequency. — X-ray brightness 
oscillations have been observed during X-ray 
bursts in four sources in which kilohertz QPO 
pairs have been detected (see ]28|jn| |). The fre- 
quency t'burst of the burst oscillation appears to 
be the stellar spin frequency I'spm or twice v S pin 
|36|]3^ ]. In these four sources, the separation Av 
between the kilohertz QPOs appears to be t'burst 
or half I'burst- 

Photon-energy dependence. — The relative am- 
plitudes of many kilohertz QPOs increase steeply 
with photon energy in the 2-15 keV energy band. 

In order to be successful, a model of the kilo- 
hertz QPOs must explain these characteristics 
and be consistent with the data obtained previ- 
ously on the 2-20 keV X-ray spectra and 1-100 Hz 
X-ray variability of the Z and atoll sources. We 
now review the proposed models, comparing them 
with these requirements. 



3. KILOHERTZ QPO MODELS 

3.1. Photon bubble oscillations 

Photon bubble oscillations (PBOs) were first 
proposed as a model for X-ray brightness varia- 
tions in accretion-powered pulsars (see ]l4| and 
references therein), where the strong stellar mag- 
netic field funnels the flow onto the magnetic po- 
lar caps, producing a mass flux onto the caps 
that is locally super-Eddington. Under these con- 
ditions, photon "bubbles" form in the accretion 
column above the caps, rise, and burst, causing 
brightness variations. In the initial simulations 
[|]~4fj , which considered an accretion rate of 0.2 Me 
and a surface dipole magnetic field of 3 x 10 12 G, 
brightness fluctuations occurred over a narrow 
range of frequencies at a few thousand Hertz. 

Klein et al. (l^] suggested that the same mech- 
anism is responsible for the kilohertz QPO pairs, 
and predicted that additional QPOs would be 
found at higher frequencies. Observations of 
Sco X-1 performed to check this jl4| showed no 
evidence of higher-frequency oscillations down to 
fairly small amplitudes. However, if there is scat- 
tering material around the neutron star magneto- 
sphere, as inferred from EXOSA T and Ginga ob- 
servations (see Fig. 1), the observed amplitudes 
of PBOs with frequencies above a kilohertz would 
be significantly reduced. In more recent simu- 
lations carried out for a magnetic field of 5 x 
10 9 G, one or two QPOs with rms amplitudes up 
to 3% occurred, with frequencies that increased 
with increasing accretion rate, accompanied by 
weaker higher-frequency QPOs [Q. The rele- 
vant three-dimensional radiation-hydrodynamics 
is very complicated, making it hard to achieve an- 
alytic insights (but see Q ) and difficult to survey 
the full parameter space numerically. 

The locally super-Eddington mass fluxes re- 
quired for photon bubbles to form are not ex- 
pected in the atoll sources, which are thought 
to have weak magnetic fields and total accretion 
rates as low as ~ 0.003 Me, but suitable condi- 
tions may exist in the Z sources. If the Z-source 
kilohertz QPOs are PBOs, then the atoll-source 
kilohertz QPOs, which appear very similar (see 
§2), would have to be generated by a different 
mechanism. It remains to be seen whether PBOs 
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can explain key characteristics of the kilohertz 
QPO pairs, such as their constant frequency sep- 
aration and the steep increase in their relative 
amplitude with increasing photon energy. 

3.2. Disk oscillations 

Titarchuk & Muslimov |39j have proposed that 
the frequencies of the burst oscillation and the 
two kilohertz QPOs seen in many sources (see 
§2) are the frequencies of three specific oscilla- 
tion modes of the accretion disk. In this model, 
as in all models in which the kilohertz QPOs are 
produced by disk oscillations, it is difficult to un- 
derstand how QPOs with amplitudes as large as 
15% of the total countrate can be generated, be- 
cause typically < 20% of the accretion luminosity 
is released in the accretion disk. In addition, one 
must explain why only three of the large number 
of disk oscillation modes produce large-amplitude 
brightness oscillations. 

No proposal was made by Titarchuk & Mus- 
limov for how the modes they discussed are ex- 
cited, why only the modes at a particular radius 
are excited, or how they produce large-amplitude 
brightness oscillations. Moreover, the expression 
for mode frequencies that they used generally pre- 
dicts that the ratio of the frequencies of two dif- 
ferent modes will remain constant as the frequen- 
cies change, whereas it is the separation between 
the kilohertz QPOs in a pair that is nearly con- 
stant. For the frequency separation to be con- 
stant in this model, one must appeal either to 
a coincidence or to some physical principle that 
has not yet been elucidated. A more serious dif- 
ficulty is that the method Q used by Titarchuk 
&; Muslimov to compute the frequencies of os- 
cillation modes in the centrifugally supported, 
strongly sheared accretion disk flow is a pertur- 
bation method designed for computing the small 
changes in the frequencies of the nonradial oscilla- 
tion modes of a spherical, gas-pressure-supported 
star caused by slow, uniform rotation. 

3.3. Magnetospheric Keplerian frequency 
and beat-frequency oscillations 

The magnetospheric beat-frequency model was 
developed to explain the single, ~15-60 Hz 

"horizontal branch oscillation" (HBO) observed 



in the Z sources (see §1). In this model, the 
frequency of the HBO is the difference between 
the Keplerian orbital frequency n m at the main 
radius where the stellar magnetic field picks up 
and channels gas from the accretion disk onto the 
magnetic polar regions (the "main gas pick-up ra- 
dius") and the stellar spin frequency Z4 p j n . In 
applying the magnetospheric beat-frequency idea 
to the kilohertz QPO pairs, Strohmayer et al. ]35[ | 
interpreted the frequency of the higher-frequency 
QPO in a pair as VKm and the frequency of the 
lower-frequency QPO as VKm — v sp i n . 

An attractive feature of the magnetospheric 
beat-frequency interpretation of the kilohertz 
QPO pairs is that it explains naturally why the 
frequency separation between the QPOs in a pair 
is nearly constant in most sources and equal to the 
burst oscillation frequency or half this frequency. 
However, there are many serious difficulties with 
this interpretation plfl . For example, no mech- 
anism has been found that will generate a rela- 
tively coherent QPO at the orbital frequency at 
the main gas pick-up radius. The high observed 
coherence of the kilohertz QPOs is also a serious 
difficulty for the magnetospheric beat-frequency 
interpretation. HBOs and kilohertz QPO pairs 
occur simultaneously in at least five of the six 
known Z sources (see Q), so the magnetospheric 
beat-frequency model cannot explain both types 
of QPO. The fact (see Fig. |) that HBOs are seen 
only in the Z sources, which are thought to have 
magnetic fields ~ 10 9 -10 10 G, and have higher 
amplitudes in the Z sources that have stronger 
magnetic fields, whereas the kilohertz QPOs are 
much stronger in the atoll sources, which are 
thought to have magnetic fields ~ 10 7 -5 x 10 9 G, 
is strong evidence that the HBOs are magneto- 
spheric and that the kilohertz QPOs are not. 

3.4. Sonic-point Keplerian frequency and 
beat-frequency oscillations 

The sonic-point model |28| was developed 
specifically to explain the kilohertz QPO pairs. 
In this model, some accreting gas spirals inward 
in nearly circular Keplerian orbits until it is close 
to the neutron star, where radiation forces or gen- 
eral relativistic effects cause a sudden increase in 
the inward radial velocity, which becomes super- 
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Figure 2. Correlation of the observed amplitudes of HBOs (left panel) and kilohcrtz QPOs (right 
panel) with the strength of the neutron star magnetic field (dashed lines) inferred from detailed radiation 
transport calculations of the X-ray spectrum [30-32]. The areas of the shaded circles are proportional 
to the measured amplitudes of the QPOs. The X's indicate that no QPO has been detected. The three 
sources in the upper left are, from top to bottom, 4U 1728-34, 4U 1820-30, and 4U 1636-536. The 
four sources in the center are 4U 1735—444, GX 3+1, GX 9+9, and GX 13+1. The four sources in the 
lower right are the Z sources Sco X-l, GX 17+2, GX 5-1, and Cyg X-2. From [31]. 



sonic within a small radial distance. The sharp 
increase in the radial velocity is usually caused by 
the drag exerted on the accreting gas by radiation 
from the star, but may instead be caused by gen- 
eral relativistic corrections to Newtonian gravity 
if the gas in the Keplerian disk flow reaches the 
innermost stable circular orbit (ISCO) without 
being significantly affected by radiation. For con- 
venience, the radius at which the radial velocity 
increases is referred to as the "sonic point" , even 
though the sonic point itself is not especially sig- 
nificant in this model. 

In the sonic-point model, the frequency of the 
higher-frequency QPO in a pair is the orbital fre- 
quency z^Ks near the sonic radius. Density fluc- 
tuations at this radius are mapped onto the stel- 
lar surface, producing an X-ray beam that rotates 
around the star with frequency vks- A distant ob- 
server therefore sees quasi-periodic X-ray spectral 
and brightness oscillations with frequency r^Ks- 

The frequency of the lower-frequency QPO in 
a pair is the difference between this orbital fre- 
quency and the stellar spin frequency. The lower- 
frequency QPO is generated by weak X-ray beams 
produced by funncling of part of the accretion 



flow near the star by the star's weak magnetic 
field. These beams rotate with the star and mod- 
ulate the radiation drag acting on the gas at the 
sonic radius, modulating the inward mass flux 
and the luminosity at the sonic-point beat fre- 
quency VBs (^Ks - ^spin Or lS Ks - 2i4 pin ). 

The sonic-point model is consistent with the ac- 
cretion rates, stellar magnetic fields, and scatter- 
ing optical depths inferred previously from EX- 
OS AT and Ginga observations of the atoll and 
Z sources and accounts for the main features of 
the kilohertz QPOs (see pq]), including their high 
and variable frequencies, their high amplitudes 
and coherences, and the common occurrence of 
kilohcrtz QPOs in pairs. Like the magneto- 
spheric beat-frequency interpretation, the sonic- 
point beat-frequency model explains naturally 
why the frequency separation between the fre- 
quencies of a kilohertz QPO pair is nearly con- 
stant and equal to the burst oscillation frequency 
or half this frequency. 

The sonic-point model also explains why the 
frequency ranges of the kilohertz QPOs are sim- 
ilar in sources with very different accretion rates 
and magnetic fields (see §4). Finally, it accounts 



Figure 3. View of the neutron star and disk along the rotation axis of the disk, which is rotating 
counterclockwise in this view, (a) The spiral trajectory followed by a single element of gas as it falls 
supersonically from the sonic radius to the stellar surface, (b) The spiral pattern of higher gas density 
formed by gas streaming inward along spiral trajectories with the shape shown in (a), from a clump 
orbiting near the sonic radius. The surface density of the disk flow is much smaller inside the sonic 
radius (lighter shaded region) than outside (darker shaded region), because of the sharp increase in the 
inward radial velocity at the sonic radius. The spiral trajectory and density pattern are from fully general 
relativistic numerical computations of the gas dynamics and radiation transport in the inner disk [28]. 



(see pH ) for the observed increase in QPO am- 
plitude with increasing photon energy and the 
observed anti-correlation between kilohertz QPO 
amplitude and stellar magnetic field strength (see 
Fig. |). 

There are several unresolved questions related 
to the sonic-point model. Perhaps the most im- 
portant is what fraction of the gas in the accre- 
tion disk does not couple strongly to the stel- 
lar magnetic field and therefore spirals inward in 
nearly circular orbits until it is close to the neu- 
tron star. This question is particularly important 
for the Z sources, some of which have inferred 
magnetic fields as high as 10 10 G. Although it has 
been realized for more than two decades (see, e.g., 
||||[l0 17 1) that instabilities will allow some frac- 
tion of the accreting gas to penetrate deep in- 
side the magnetosphere, it has not been possible 
to calculate this fraction reliably or to determine 
how much of this gas remains in nearly circular 



orbits. The kilohertz QPOs are strong evidence 
that at least some gas is orbiting close to the star, 
because the orbital frequencies at the radii where 
gas is expected to begin coupling to the stellar 
magnetic field are ~5-10 times lower than the 
frequencies of the kilohertz QPOs. 

The sonic-point model is currently the most 
fully developed and successful explanation for the 
kilohertz QPOs. In the next section we describe 
some of the calculations that have been carried 
out to explore it. 

4. THE SONIC-POINT MODEL 

Figure || shows how the QPO at the sonic-point 
Keplerian frequency is generated in the sonic- 
point model. Gas in the accretion disk drifts 
slowly inward until it becomes exposed to radi- 
ation from the stellar surface. Once exposed, the 
gas loses its angular momentum to radiation drag 
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Figure 4. Results of fully general relativistic numerical computations [|28j of the gas dynamics and 
radiation transport in the inner disk in the sonic-point model, for a neutron star of radius 5M. (a) The 
inward radial velocity v r of the gas in the disk measured by a local static observer, as a function of the 
Boyer-Lindquist radial coordinate expressed in units of the stellar mass M. (b) The radial optical depth 
r r from the stellar surface through the disk flow to the radius shown on the horizontal axis. The four 
curves in each panel are labeled with the assumed accretion rate M measured in units of the accretion 
rate Me that would produce an accretion luminosity at infinity equal to the Eddington critical luminosity. 



in a radial distance Ar <0.01r and falls inward 
supersonically along spiral trajectories like that 
shown in Figure ^|a. Gas falling inward from the 
sonic radius along spiral trajectories collides with 
the neutron star around its equator, producing an 
X-ray emitting equatorial ring, which is indicated 
by the grey ring around the star. 

Magnetoturbulence, differential cooling, and 
radiation forces create density inhomogeneities 
("clumps") in the gas in the accretion disk. Gas 
streaming inward from clumps orbiting near the 
sonic radius along trajectories with the shape 
shown in Figure |3|a generates the density pattern 
shown in Figure [|b. Collision with the stellar 
surface of the denser gas from the clumps cre- 
ates beams of brighter X-ray emission, like the 
beam indicated by the white dashed lines in Fig- 
ure ||b, which move around the star's equator 
with a frequency equal to the orbital frequency 
of the clumps, generating a quasi-periodic bright- 
ness oscillation with frequency isk s . 



The QPO at the sonic-point beat frequency is 
generated by weak X-ray beams that rotate with 
the star and modulate the radiation drag acting 
on the gas at the sonic radius once each beat 
period. As a result, the inward mass flux from 
each clump oscillates, causing the luminosity of 
the footprints to oscillate at the sonic-point beat 
frequency v^ s - 

The nature of the transition to supersonic ra- 
dial inflow is illustrated by the radial velocity and 
optical depth profiles shown in Figure |J. These 
profiles were calculated assuming that the az- 
imuthal velocity of the gas in the disk is nearly 
Keplerian far from the star. Internal shear stress 
in the disk was assumed to create a constant in- 
ward radial velocity if ', as measured in the local 
static frame, of 10 -5 . The half-height h(r) of the 
disk flow at radius r was assumed to be er at all 
radii, where e is a constant and r is the radius, 
and the kinetic energy of the gas that collides 
with the surface of the star was assumed to be 
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converted to radiation and emitted from a band 
around the star's equator with a half-height equal 
to eR. For simplicity, and to show the effects 
of radiation forces more clearly, any effect of the 
stellar magnetic field on the gas dynamics near 
the sonic transition was neglected. 

Once the drag force exerted by the radiation 
from the stellar surface begins to remove angular 
momentum from the gas in the Keplerian disk, 
centrifugal support is lost and the gas falls in- 
ward, accelerating rapidly. Radiation that comes 
from near the star and is scattered by the gas in 
the disk is usually scattered out of the disk plane 
and hence does not interact further with the gas 
in the disk. Moreover, second and successive scat- 
terings do not contribute proportionately to the 
azimuthal radiation drag force on the gas because 
the radiation field is aberrated by the first scat- 
tering and afterward carries angular momentum 
]25| , p6[ | . The interaction of the radiation with the 
gas in the disk was therefore treated by assuming 
that the intensity of the radiation coming from 
the star is attenuated as it passes through the 
gas in the disk, diminishing as exp(— r r ), where 
r r (r) is the Thomson scattering optical depth ra- 
dially outward from the stellar surface to radius r, 
and that scattered radiation does not contribute 
to removal of angular momentum from the gas in 
the inner disk. In order to simply calculation of 
the radiation drag force, the differential scatter- 
ing cross section was assumed to be isotropic in 
the frame comoving with the accreting gas (see 
pot). The radiation field and the motion of the 
gas were computed in full general relativity. 

Figure || shows that in this model, the sonic ra- 
dius decreases with increasing accretion rate un- 
til it reaches the radius of the ISCO, after which 
it stops decreasing. Hence, the frequency of the 
sonic-point Keplerian frequency QPO increases 
steeply with increasing luminosity until the sonic 
radius reaches the radius of the ISCO, after which 
it stops increasing, as shown in Figure [s]. 

The sonic point must occur between the radius 
R ms of the innermost stable circular orbit and 
the radius where the two curves shown in Fig- 
ure ^| intersect, which is ~ 3i? ms . The frequency 
of the sonic-point Keplerian frequency QPO is 
therefore restricted to a relatively small interval 
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profiles shown in Figure f|. From [28] . 



for all sources. For example, if the neutron stars 
in LMXBs all have masses ~1.7M , then the fre- 
quencies of their sonic-point Keplerian frequency 
QPOs would all be confined to the range ~400- 
1300 Hz. 

5. DISCUSSION 

The discovery with RXTE of -1000 Hz bright- 
ness oscillations from a large number of accreting 
neutron stars is a spectacular achievement that 
validates both the scientific expectations that led 
to the mission and the long years of hard work 
that were needed to bring it to fruition. 

As explained in §3, the observations made to 
date strongly favor beat-frequency models of the 
kilohertz QPOs detected in the persistent emis- 
sion. The magnetospheric beat-frequency inter- 
pretation of the kilohertz QPO pairs can account 
for the fact that their separation frequencies are 
approximately constant and related to the fre- 
quencies of the brightness oscillations seen during 
X-ray bursts, but this interpretation suffers from 
many serious difficulties. The sonic-point beat- 
frequency model also explains this fact and is the 
most fully developed and successful model, but it 
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is not yet confirmed. Therefore, it is important 
to consider which further observations and cal- 
culations would be particularly helpful in testing 
the beat-frequency hypothesis and discriminating 
between these two beat-frequency models. 

The detection in the persistent emission of sta- 
ble oscillations with frequencies harmonically re- 
lated to those seen during X-ray bursts would 
confirm the spin-frequency interpretation of the 
latter and strongly support the beat-frequency in- 
terpretation of the kilohertz QPO pairs. On the 
other hand, detection of stable oscillations with 
frequencies not harmonically related to those seen 
during X-ray bursts would seriously undermine 
both models. 

Beat-frequency models also predict that weak 
oscillations will be present at other special fre- 
quencies, such as the first overtone of the beat fre- 
quency and the sum of the spin frequency and the 
orbital frequency, and detection of any of these 
would provide strong support for these models. 
It is therefore important to search for QPOs at 
these frequencies. Even if none are detected, up- 
per limits on their amplitudes would be valuable, 
because they would constrain models of the coro- 
nae around these neutron stars (see Fig. 0). 

The sonic-point model predicts that the 
stronger the stellar magnetic field, the weaker 
the kilohertz QPOs will be (28). In the magneto- 
spheric interpretation, one expects just the oppo- 
site. Figure || shows that the current data on kilo- 
hertz QPO amplitudes clearly favors the sonic- 
point model, but more precise and uniform am- 
plitude measurements, as well as further progress 
in modeling the X-ray spectra of the Z and atoll 
sources, would help to clarify the situation. 

Besides their intrinsic interest, the ~1000 Hz 
oscillations discovered with RXTE provide a new 
tool with which to probe the nature of strong 
gravitational fields and the properties of dense 
matter. Modeling of general relativistic effects 
on the gas dynamics and radiation transport 
processes involved in the generation of kilohertz 
QPOs may provide evidence for the existence of 
innermost stable circular orbits J28|,|2§] and frame- 
dragging p5| , p6| , ^3|j34|l , both of which are impor- 
tant predictions of strong-field general relativity. 

If the frequencies of the highest-frequency kilo- 
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Figure 6. Comparison of r/f{ ow , the fraction of the 
specific angular momentum of an element of gas 
in Keplerian circular orbit at Boyer-Lindquist ra- 
dial coordinate r (measured in units of the stellar 
mass) that must be removed in order for the gas 
to fall from r to the radius R ms of the innermost 
stable circular orbit, with ?7 ra d(max) the largest 
fraction of the specific angular momentum of an 
element of gas that can be removed by radiation 
coming from the surface of a nonrotating, isotrop- 
ically radiating, spherical star of radius 5M. 



hertz QPOs are orbital frequencies, as in beat- 
frequency models, these QPOs provide impor- 
tant new constraints on the masses and radii of 
the neutron stars in LMXBs and on the equa- 
tion of state of neutron star matter [ p^|j2^j29| , ^5[ | . 
The brightness oscillations observed during X-ray 
bursts may constrain the compactness of neutron 

stars HjMi. 

New RXTE observations are continuing to 
yield significant fresh insights, and the rapid pace 
of important new discoveries is therefore likely to 
continue for many years. 
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grant NAG 5-2925 at the University of Illinois, 
and NASA grant NAG 5-2868 at the University 
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